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^ ■ The polarization of the Cosmic Microwave Background (CMB)is a powerful ob- 

■ servational tool at hand for modern cosmology. It allows to break the degeneracy 

Q_|. of fundamental cosmological parameters one cannot obtain using only anisotropy 

Q i data and provides new insight into conditions existing in the very early Universe. 

^ I Many experiments are now in progress whose aim is detecting anisotropy and 

^ ■ polarization of the CMB. Measurements of the CMB polarization are however 

^ ■ hampered by the presence of polarized foregrounds, above all the synchrotron 

"-j . emission of our Galaxy, whose importance increases as frequency decreases and 

I dominates the polarized diffuse radiation at frequencies below ~ 50 GHz. In the 

I past the separation of CMB and synchrotron was made combining observations 

of the same area of sky made at different frequencies. In this paper we show that 
the statistical properties of the polarized components of the synchrotron and dust 
foregrounds are different from the statistical properties of the polarized compo- 
nent of the CMB, therefore one can build a statistical estimator which allows 
to extract the polarized component of the CMB from single frequency data also 
when the polarized CMB signal is just a fraction of the total polarized signal. 
This estimator improves the signal/noise ratio for the polarized component of the 
CMB and reduces from ~50 GHz to ^20 GHz the frequency above which the 
polarized component of the CMB can be extracted from single frequency maps 
of the diffuse radiation. 
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1 Introduction 



This year is a decade since the first detection of the anisotropy of the Cosmic 
Microwave Background at large angular scales (> 10°) [1], [2]. Today the CMB 
anisotropy (CMBA) has been detected also at intermediate (~ (1° — 10°)) and 
small angular scales (< 1°), so the CMBA angular spectrum is now reasonably 
known down to the region of the first and second Doppler peaks [3], [4], [5]. 
Its shape gives information e.g. on the spectrum of the primordial cosmological 
perturbations or can be used to test the inflation theory but rises new questions to 
which CMBA cannot answers. Responses can on the contrary be obtained looking 
at the CMB polarization (CMBP) produced by Thomson scattering of CMB 
photons on the matter anisotropics at the recombination epoch. In particular 
one can hope to use CMBP to disentangle the effects of fundamental cosmological 
parameters like density of matter, density of dark energy etc., effects anisotropy 
do not separate.This is among the goals of space and ground based experiments 
like [6], [7], [8], [9], [10], [11], [12] and is the main goal of SPOrt a polarization 
dedicated ASI/ESA space mission on the International Space Station [13]. The 
relevance of the CMB polarization was remarked for the first time by M. Rees 
[14]. Since him many models of the expected features of the CMBP have been 
pubhshed (see for instance [15], [16], [17]). They stimulated the search for CMBP, 
but in spite of many attempts so far no CMB polarization has been detected. 
This observation is in fact extremely difficult because the expected signal is at 
least an order of magnitude smaller than the amplitude of the CMBA. Moreover 
foregrounds and their inhomogeneities cover the polarized fraction of the CMB or 
mimic CMBP spots, making the signal to noise (CMBP / polarized foreground) 
ratio unfavorable. In this paper wc will concentrate on methods for improving 
this ratio and for disentangling CMBP and polarized foregrounds. 
In the microwave range the galactic foregrounds include: 

• synchrotron radiation (strongly polarized), 

• free- free emission (polarization negligible), 

• dust radiation. 

Because here we are interested in polarization, in the following we will neglect 
the free-free emission, whose expected level of polarization is negligible. More- 
over the polarized signals produced by dust, if present, (e.g. [18], [19]), may be 
treated as an addition to the synchrotron effects. In fact, as it will appears in the 
following, the important quantities in our analysis are the statistical properties of 
the foreground spatial distribution and, by good fortune, the spatial distribution 
of the dust polarized emission is similar to that of the synchrotron emission, since 
behind both radiation types there is the same driving force, the galactic magnetic 
field which alignes dust grains and guides radiating electrons. 
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Separation of foregrounds and CMBA was successfully solved when the CMB 
anisotropy was discovered [2]. When we go from CMBA to CMBP the separa- 
tion of foreground and background however is more demanding and the problem 
of discriminating foreground inhomogneities from true CMB spots severe. Ap- 
proaches used in the past e.g. [20], [21], [22], [23], [24]) were essentially based on 
the differences between the frequency spectra of foregrounds and CMB, therefore 
require multifrequency observations. 

In this paper we suggest a different method which takes advantage of the fact 
that the measured values of the parameters we use to describes the polarization 
of the diffuse radiation when measured at a given frequency in different directions 
behave as stochastic variables. Because the mathematical and statistical proper- 
ties of these variables for synchrotron and CMB are different, we suggest to use 
statistical methods for analyzing single frequency maps of the diffuse radiation 
and disentangling their main components, synchrotron and CMBP. This method 
was proposed and briefly discussed in [25]. Here we present a more complete 
analysis. 



2 Polarization Parameters 
2.1 The Stokes Parameters 

Convenient quantities commonly used to describe the polarization status of radi- 
ation are the Stokes parameters (see for instance [26], [27], [28]). Let's assume a 
monochromatic, plane, wave of intensity I and amplitude oc \fl. In the observer 
plane, orthogonal to the direction of propagation of the electromagnetic wave 
emitted by the electron, we can choose a pair of orthogonal axes / and r. On 
that plane the amplitude vector of an unpolarized wave moves in a random way. 
On the contrary it describes a figure, the polarization ellipse, when the wave is 
polarized. 

Projecting the wave amplitude on / and r we get two orthogonal, linearly 
polarized, waves of intensity Ii and Ir,{I = Ii + Ir) whose amplitudes are oc \/Ti 
and oc y/T^ respectively. If the original wave of intensity / is polarized, // and Ir 
are correlated: let's call I\2 and /21 their correlation products. 

By definition the Stokes parameters are the four quantities: / = Ii + Ir, 
Q = Ii- Ir, U = I12 + /21, and V = i^hi - /12). 

Q and U describe the linear polarization, V the circular polarization and / 
the total intensity. 
The ratio 

tan2x=^ (1) 
gives the angle x between the vector / and the main axis of the polarization ellipse 
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(0 < X < tt). 

Rotating the {I, r) coordinate system by an angle we get a new coordinate 
system {l,f) in which the Stokes parameters become 

Q = Q cos 2(j) — U sin 20 

(j = [/ cos 20 + (5 sin 20 

When = X = (1/2) arctan^ the axes of the polarization eUipse coincide with 
the reference axes / and r (see fig.l). 

2.2 Electric and Magnetic Modes 

To analyze the properties of the CMB polarization it is sometimes convenient 
to use rotationally invariant quantities, like the radiation intensity / and two 
combinations of U and Q: Q+iU and Q — iU. The intensity / can be decomposed 
into usual (scalar) spherical harmonics Yi^{9,(p). 

I = ^airaYlm{0-,(p) (3) 

The quantities Q ±iU can be decomposed into ±2 spin harmonics [29], [30], [31] 

g±iC/ = ^a,t2Y,±^(^,(^) (4) 

l,m 

The ±2 spin harmonics form a complete orthonormal system (see, for instance, 
[32], [33], [34], [35]) and can be written [30], [36]: 

(5) 

r,-2(^,(^) = 7V-2p^-2(^)e-v= 



^Alternatively, one can use the equivalent polynomials derived in [32] 
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where 

pLi^) = (1 - + x)-^p^rr''-""{x) (6) 

is a generalized Jacobi polynomial, s = ±2 and: 



(m + s) is — m) , , ^ 



{l-s)\{l + s)\ 



'"^ 2^V 47r ^ (Z-m)!(Z + m)! 

is a normalization factor. 

The harmonics amplitudes af^ correspond to the Fourier spectrum of the an- 
gular decomposition of rotationally invariant combinations of Stokes parameters. 

Because spin ±2 spherical functions form a complete orthonormal system: 

/ ^)Y;^^{e, ip)dQ = SimSvm' (7) 

we can write 

af^ = f dn {Q{e, cp) ± iu{0, cp)) y;^^' (s) 

Following [31] we now introduce the so called E (electric) and B (magnetic) 
modes of these harmonic quantities: 

(9) 

They have different parities. In fact when we transform the coordinate system 
Oxyz into a new coordinate system Oxyz, such that 



/ = / 

— * — * 

r — —r 

the E and B modes transform in a similar way: 

d^ = 

Q remains identical in both reference systems and U changes sign. 
It is important to remark that and are uncorrelated. 
In terms of Q and U we can write: 

at^y dn (O (Y,i' + 1^-^) + iU {Y,t' - 1^-^)) 

ag^^^jdn (iQ(y+2-i^-^)-c/(y+^ + i^-^)) 



(10) 



(11) 



(12) 
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therefore : 

{(aLr) - {{at?) = 2jdn {{m) - {\U'\)) {Y,l'YC-' + YC^^-J) (13) 

Here and designate values of correlators of delta correlated 2D 

stochastic fields Q and U. Omitting mathematical details, the correlation equa- 
tions for Q and U are: 

{QQ*) = \Q^\S{n-n') 

(UU*) = \U''\S{Q-Q') (14) 
{QU*) = 

where S{Q — Q') = S{cos6 — cos6') ■ S{ip — </?') is the Dirac delta - function on the 
sphere. (In the following we will sometimes omit indexes I and m). 



3 Synchrotron Radiation and its Polarization 

Synchrotron radiation results from the helical motion of extremely relativistic 
electrons around the field hues of the galactic magnetic field (see, for instance 
[26], [28], [37]). The electron angular velocity 

eHp rrieC^ rUeC^ 

OJe = ^ = t^o— ^ (15) 

nieC t t 

is determined by the ratio between Hp, the component of the magnetic field 
orthogonal to the particle velocity, and £, the electron energy. As it moves 
around the magnetic field lines the electron radiates. 

a) Single electron 

Until the circular velocity is small {v « c, cyclotron radiation) the electron 
behaves as a rigid dipole which rotates with gyrofrequency (15) in a plane or- 
thogonal to the magnetic field direction and emits a single line. The spatial 
distribution of the radiation has dumbbell shape (see fig. 2): 

J(e,$) ~ (1 + cos^0) (16) 

The radiation is circularly polarized along the dumbbell axis (G = 0) and linearly 
polarized in directions orthogonal to it (0 = 90°). 

When the electron velocity increases the radiation field changes until at f c, 
{S ^ mgC^) it assumes the peculiar charachters of synchrotron radiation: 

i) radiated power proportional to f ^ and H^, 
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Figure 2: Spatial distribution of the cyclotron radiation produced by a single electron 




Figure 3: Velocity cone of an ultrarelativistic {y ~ c) electron {H = magnetic vector, 
k wave vector, V angle between the electron velocity and the direction of observation) 

ii)continous frequency spectrum so concentrated around: 



<^m = t^o OC Hp£^ (17) 

1 cos 1p 

c 

(■0 is the angle between the velocity vector v and the wave vector k, see fig. 3), 
that in a given direction emission can be assumed monochromatic, 

iii)radiation almost totally emitted in the forward direction of the electron 
motion, inside a narrow cone ^ of aperture (see fig.4) 

^ ~ (18) 
Inside that cone (cos ip ^ 1) the frequency is maximum and equal to 

\m.eC J 



^the symmetry plane(0 = 90°) of the cyclotron dumbbell beam, seen by a fast moving 
observer [v « c) becomes a cone folded around the direction of movement 
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Figure 4: Distribution of the radiation emitted by an ultrarelativistic {v ~ c) electron 
with istantaneous velocity v spiralling (dashed line) around the lines of the magnetic 
field H (orthogonal to the sheet) 

In the opposite direction (cosV' ~ —1) intensity and frequency are sharply re- 
duced. 

iv)radiation 100 % linearly polarized when seen along the surface of the emis- 
sion cone.Inside the cone the linear polarization is still dominant but a small 
fraction of circular polarization exists, {V ~ ci(^^^^) and I ^ U Q). Outside 
the cone the very small fraction of radiation produced is elliptically polarized and 
becomes circularly polarized when seen in directions orthogonal to the circular 
component of the electron spiral motion, i.e. along the Hp direction. So the 
Stokes parameters depend on H, the angle // between the H and the hue of sight, 
and the dimensionlcss frequency u/uc, where Uc = 1.5{u!m,o/^'^) is the so called 
critical frequency, [26], [28], itself function of H (see eq.l5). 

b) Cloud of monoenergetic electrons 

When the effects of many monoenergetic electrons with uniform distribution of 
pitch angles are combined, I, Q and U are reinforced (the Stokes parameters are 
additive) while V is erased. In fact 

= c,HpF(fJ 
Q{u) ^ C2HpFp{f) cos2x 

(20) 

U{iy)=C2HpFp{fJsm2x 
y(i/) fti 

where ci and C2 are constants, F and Fp frequency functions (almost monochro- 
matic), and X the angle between the projection of the magnetic vector on the 
observer plane and an axis on that plane (The projection of the magnetic vector 
on the observer plane is the minor axis of the polarization ellipse). When the 
magnetic field direction varies, the angle x varies therefore cos 2% and sin 2% must 
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be averaged on the magnetic field distribution along the fine of sight through the 
cloud of emitting electrons. In conclusion the degree of polarization p varies 
between a maximum value, (uniform magnetic field) and zero (magnetic field 
randomly distributed). 

c) Electrons with power law energy spectrum 

In the interstellar medium the radiating electrons are the cosmic ray electrons 
whose spectrum is a power law energy spectrum (see for instance [38], [39]) and 
references therein): 

N{S) = )C£-^ (21) 

where 7 ~ 2.4 — 3.0 is the spectral index, and K a normalization constant. 
Because the radiation produced by each electron is practically monochromatic 
the resulting radiation spectrum is also a power law: 

= Io{^)H^u-^' (22) 

where P' = ^ is the intensity spectral index, P — P'+2 the temperature spectral 
index and /q is a slow function of 7 [26]. When the magnetic field is not uniform 
Hp is replaced by < Hp > and Jq by a slightly different function of 7. 
The Stokes parameters, products of the intensity /(t^), the degree of polarization 
p and cos 2% or sin 2%, are: 

7+1 _2-l 

Q ^ cs Hp'' u 2 cos 2x 

[/ = C3 H^ u-"^ sin 2x 
V ^0 

Here C3 is a constant determined by 7 and by the magnetic field distribution 
along the integration path. Faraday rotation, induced by the combinations of 
thermal electrons (if present) and magnetic field along the line of sight through 
the synchrotron emitting region, usually reduces the polarization level. Moreover 
when the field is not uniform we have to use < sin 2% > and < cos 2% > instead 
of sin 2% and cos 2% so the degree p of polarization decreases. Summarizing p 
varies between (random magnetic field distribution) and: 

Pmax = < 1 (24) 

07 + / 

(uniform magnetic field and Faraday rotation absent). 

When one looks in different directions through the interstellar medium H, 
Hp — H sin fjL, X, 7, P and K vary (see for instance [40], [41]. Practically their 
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values, averaged along the line of sight, behave as random variables ^ therefore 

the Stokes parameters Qg and Ug (s stays for synchrotron) associated to the 
galactic (synchrotron) foreground behave as stochastic functions of the direction 
of observation. We can therefore write: 

(Q^) = {Ul) (25) 

It follows that for synchrotron radiation if 7^ also Us ^ 0, therefore both 
electric and magnetic modes exist: 

a^'' 7^ 0, a^'^ 7^ (26) 

moreover for the synchrotron background we can write (see eqs. (13) and (25)): 

{{a^'r) = {{a^n') (27) 

4 CMB and its Stokes Parameters 

In a homogeneous and isotropic Universe the only quantities which change as the 
Universe expand are temperature and intensity I — Ii + Ir'- both decrease adia- 
batically. Because this is true for Ii and separately, we do not expect anisotropy 
nor polarization therefore Q = and U = /„ = is a natural consequence. 

On the contrary, inhomogeneities and perturbations of the matter density or 
of the gravitational field, induce anisotropy and polarization of the CMB. At the 
recombination epoch linear polarization is produced by the Thomson scattering of 
the CMB on the free electrons of the primordial plasma. The polarizarion tensor 
can be calculated solving the Boltzman equations, which describe the transfer of 
radiation in a nonstationary plasma permeated by a variable and inhomogeneous 
gravitational field [42], [43], [44], [15]. 

In our Universe the gravitational field can be divided in two parts: a back- 
ground field, with homogeneous and isotropic FRW metric, and an inhomoge- 
neous and variable mix of waves: density fluctuations, velocity fluctuations, and 
gravitational waves. Because of their transformation laws these waves are also 
said scalar, vector and tensor perturbations, respectively. 

Scalar (density) perturbations affect the gravitational field, the density of 
matter and its velocity distribution. They were discovered by astronomers who 
studied the matter distribution in our Universe on scale from ~ 1 Mpc to ~ 100 
Mpc. It is firmly believed they are the seeds of the large scale structure of the 
Universe and are reflected by the large scale CMB anisotropy detected for the 
first time at the beginning of the '90s [1], [2]. Their existence is predicted by the 
great majority of models of the early Universe. 

^e.g. F = [H sin jj) ^ is a nonzero mean, nonzero variance variable; cos 2% and sin 2% are 
zero mean, nonzero variance variables 
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Vector perturbations, associated to rotational effects, perturb only velocity 
and gravitational field. They arc not predicted by the inflation theory and it is 
common believe that they do not contribute to the anisotropy and polarization 
of the CMB. 

Tensor perturbations affect exclusively the gravitational field. 

In the following we will concentrate our attention on density waves, the dom- 
inant perturbation. Becaiisc observation shows that the effects of these pertur- 
bations are small, we can treat them as small variations 5/, 5^., 5„ of and /„. 
If we introduce the auxiliary functions a and 13: 

5i-5r = {l- /x^)/?, (28) 
<^« = 0, 

{fi is the angle between the line of sight and the wave vector) and assume plane 
waves, the Boltzman equations (see for instance [15], [29], [30] and reference 
therein) become: 

^ = F - ^aTnea{r])a - ^aTnea{r])p 

(29) 

^ ^ -ji](^Tnea{r))a - ^aTnea{r))(3 

where F is the gravitational force which drives both anisotropy and polarization, 
(7t is the Thomson cross-section, Ug is the density of free electrons, and 0(77) the 
scale factor. 

These equations give: 

Q = -^(1 - /x^) / F{n) (e-- - e-^-) drj 

(30) 

C/ = 

where T{'r]) is the optical depth of the region where the phenomenon occurs. 
Rotating the coordinate system we can generate a new pair of Stokes parameters 
{Q', U'): no matter which is the system of reference we choose these parameters 
satisfy the symmetry parity condition. 

Because as we saw for the CMB there is a system in which Q 7^ and U — 0, 
we may conclude that magnetic modes of the CMB polarization vanish and only 
electric modes exist [31]. Therefore for the CMB we can write: 

where index d stays for density perturbation. 
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Concluding: in our approximation the CMB polarization is a random tensor 
field whose characteristcs are set by primordial density perturbations (see, for 
instance [45]). 

5 Separation of the polarized components of Syn- 
chrotron and CMB Radiation 

The CMB radiation we receive is a diffuse background which reaches us mixed 
with foregrounds of local origin. When the anisotropy of the CMB was detected, 
to remove the foregrounds from maps of the diffuse radiation, data were reorga- 
nized in the following way: 

Td{x, y) = fn{x, y) + Mi{x, y)fi^c{x, y) (32) 

where T^. is a two dimension vector (map) which gives the total signal measured 
at different points {x^y) on the sky, T„ the noise vector, M the matrix which 
combines the components Tj,c of the signal. At each point (x, y) we can in fact 
write: 

-^d 9n^n ~l~ 9cmb^cmb,c ~l~ 9syn^syn,c ~l~ 9ff'^ff,c ~l~ 9dust^dust,c (33) 

where Qi arc weights, given by M, Tsyn,c is the synchrotron component, Tjjc the 
free- free emission component, Tdust.c the dust contribution and so on. Using just 
one map the signal components cannot be disentangled. If however one has maps 
of the same region of sky made at different frequencies it is possible to write a 
system of equations. Provided the number of maps and equations is sufficient, 
the system can be solved and the components of separated, breaking the 
degeneracy. We end up with a map of T^^h which can be used to estimate the 
CMB anisotropy. 

When we look for polarization at each point on the sky we measure tensors 
components instead of scalar quantities, therefore to disentangle the polarized 
components of the CMB we need a greater number of equations. Here we will 
concentrate on the separation of the two dominant components of the polarized 
diffuse radiation: galactic synchrotron (plus dust) foreground and CMBP. 

5.1 The estimator D 

Instead of observing the same region of sky at many frequencies, we suggest a 
different approach. It takes advantage of the differences between the statistical 
properties of the two most important components of the polarized diffuse radia- 
tion (CMB (background) and synchrotron (foreground) radiation) and does not 
require multifrequency maps. 
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We define the estimator: 

D = {{a^f) - {{ay) (34) 

where — a^'^ + a^''^ and — a^'^ + a^''^ (here and in the following indexes s 

or d stay for synchrotron and density perturbations, respectively). 

Because 

• E and B modes of synchrotron do not correlate each other neither correlate 
with the CMB modes 



(a^)2) = ((a^'^)2) + ((a^''^)2) + 2 (a^'^a^-'^) = ((a^'^)^) + ((a^''^)^), 

• Electric and Magnetic modes of the CMB polarized component satisfy the 
condition 

a^''^ ^ a^^'^ = 

• Electric and Magnetic modes a^'^ and a^'^ . of the synchrotron radiation 
are both different from zero and their average values identical, 



we get 



D = ((a^-'^)') (35) 

It means that D provides an estimate of the contribution of the CMB to maps of 
the polarized diffuse radiation contaminated by galactic synchrotron emission. 

Let' s now consider the angular power spectrum of D. For multipole I we can 
write: 

A = {afr - {aff = ^ E {\aL? ' \aL?) (36) 

m=—l 

where laf^]"^ and |a^P are random variables with gaussian distribution p{af^) 
(see eqs.(54) and (55) in Appendix A). According the ergodic theorem (in the 
limit of infinite maps, the average over 2D space is equivalent to the average 
over realisations) the average value of Di is equal to the difference of the average 
values of |a;f„P and |a^|^ summed over m. Taking into account equation (55) we 
can therefore write: 

(A) = (|afp-|afr) (37) 

where ^: 

''Equation (38) is an explicit form of the average of the stochastic variables and |a^P 

over a probability density piaf^^), the short form being triangle brackets 
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—oo 

(38) 

oo , , 

(aff = / lafMafJd'^at = (\at?) 

— oo 

Comparing eq.38 with the ordinary definition of multipole coefficients: 

Cr = ^, E (l«f/l') (39) 

m=—l 

it immediately appear that we can write: 

(40) 

(afy = CP 

5.2 Separation uncertainty 

For synchrotron radiation {Di) should be zero, non zero for density perturbations, 
but on real maps it is always different from zero. In fact a map is just a realization 
of a stochastic process and the amplitudes of |a^p and |a^„p, averaged over 2D 
sphere, have uncertainties which add quadratically making (Di) ^ even in the 
case of synchrotron polarization. This effect, very similar to the well known 
"cosmic variance" of anisotropy [46], [47] (the real Universe is just a realisation 
of a stochastic process, therefore there will be always a difference between the 
realization we measure and the expectation value) does not vanish if observations 
are repeated. 

The variance of Di is 

V{Df) = (Df) - {D,f (41) 

Being sums over m, {—l<m< +/) of 2/ + 1 stochastic values with gaussian 
distribution, and have distributions with 2 * (2/ + 1) degrees of freedom, 
therefore we expect their variances can be written as 

^(«"'')^ « («' 

More explicitly 

(A^) = ^((«fr+(«fr) (43) 

and when the synchrotron foreground is dominant (of = of) 

- ^(«f (44) 

in agreement with (42). 
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5.3 A criterium for CMBP detection 



The synchrotron foreground is a sort of system noise which hampers detection 
of the signal, the CMB polarization. At frequencies sufficiently high (fs above 
50 GHz, as we will sec in the next section) the noise is small compared to the 
signal therefore direct detection of CMBP is possible. At low frequencies on the 
contrary the CMBP signal is buried in the synchrotron plus dust noise. In this 
case to recognize the presence of the CMB polarization we can use our estimator 
D. 

At angular scale I, to be detectable the CMBP must satisfy the condition 

Cf ''^ > A ■ Cf'' (45) 

where C^^'* are the coefficients of the multipole expansion of the E modes and A 
is the confindence level of the signal detectability. In a similar way we can write 
for our estimator: 

Df''^ > A ■ Df'' (46) 
where, as we said before, Df'^ = Cf'^ and 

Df^' = /^(«f (47) 
Therefore, the criterium for the CMBP detection becomes 

Df^' > A ■ ][^Cf^' (48) 

This means that using D we can recognize CMBP in a map of the polarized 
diffuse radiation with an uncertainty ~ \j 2l \ l ^^'^ which decreases as I and 
angular resolution increases. 



6 Angular power spectra of polarized synchrotron. 

The angular power spectra of the polarized component of the synchrotron radi- 
ation have been discussed by [48], [49], [50], [51] and [52]. It appears that the 
power spectra of the degree of polarization p and of electric and magnetic modes 
E and B follow power laws up to Z ~ 10^. The spectral index of the power law 
which holds for the degree of polarization is ~ 1.6 — 1.8. For the E and B 
modes different authors get from the observations values of the spectral index 
different, but marginally consistent. In paper [51] and [52], using Parkes data, 
the authors get 
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with ue ~ 0(B ~ 1.4-=- 1.5 and dependence of a on the sky region and the 
frequency. 

In paper [48], using Effelsberg and Parkes data, the authors get: 



with a — (1.8 ± 0.3), (3 — 2.9 and Cq = (1.6 ± 1) (here we modified the original 
expression given in [48] writing it in adimensional form). As said above the 
spectral indexes a obtained by [51] [52] and [48] are marginally consistent. 

Both authors got their results analyzing low frequency data (1.4, 2.4 and 2.7 
GHz), therefore the extension of their spectra to tens of GHz, the region where 
CMB observations are usually made, depend on the accuracy of (3, the temper- 
ature spectral index of the galactic synchrotron radiation. A common choice is 
(3 — 2.9 but in literature there are values of /? ranging between 2.5 and ~ 3.5. 
Moreover (3 depends on the frequency and the region of sky where measurements 
are made (see [38], [53], [54], [55]). Finally the values of P in literature have been 
obtained measuring the total (polarized plus unpolarized) galactic emission. In 
absence of Faraday effect (3, /Spoi and /Sunpoi, the spectral indexes of the total, po- 
larized and unpolarized components of the galactic emission, are identical (see 
eqs. (22), (23)). However when Faraday effect with its z/^^ frequency dependence 
is present, we expect that, as frequency increases, the measured value of the 
degree of polarization (se eq. (24)) increases, up to 



therefore we should measure f3poi < f3. The expected differences are however well 
inside the error bars of the available data therefore at present we can neglect 
differences and set (3 ~ (3poi — (3unpoi- 

Instead of extrapolating low frequency results it would be better to look for 
direct observations of the galactic emission and its polarized component at higher 
frequencies. Unfortunately above 5 GHz observations of the galactic synchrotron 
spectrum and its distribution are rare and incomplete. At 33 GHz observations 
by [56] give for some patches on the sky a galactic temperature of about 2 -^4/^i^' 
from which follows that at the same frequency we can expect polarized foreground 
signals up to several /iK. At 14.5 GHz observations made at OVRO [57] give 
synchrotron signals of 175 fiK, equivalent to about 15 /iK at 33 GHz of which 
up to 10 fiK can be polarized. 

In conclusion there are large uncertainties when one try to decide at which 
frequency the polarized component of the galactic diffuse emission does not dis- 
turb observation of the GMB. To be on the safe side we can say, and all the 
authors of papers [48], [49], [50], [51] and [52] agree, that the CMB polarized 
component surely overcomes the polarized synchrotron foreground only above 50 
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Figure 5: Power spectrum versus multipole order I of the polarized components of 

CMB (solid line) and three possible power spectra of the galactic synchrotron radiation 
(a = 1.5 (dotted line), a = 1.8 (dashed line) and a = 2.1 (dot-dash line)) calculated 
at 43 GHz. (See text for details of calculations and model) 

GHz, therefore direct observations of the CMB are better made at frequencies 
greater of 50 GHz. Below 50 GHz observations of CMBP are problematic because 
the contamination by galactic polarized emission can be high and its evaluation, 
usually by multifrequency observations, not sure. 

This conclusion is supported by figure 5. It shows plots versus the multipole 
order / of the power spectra of the polarized components of CMB and galactic 
synchrotron calculated at 43 GHz. The CMB spectrum has been obtained using 
CMBFAST [58] and standard cosmological conditions (CMB power spectrum 
normalized to the COBE data at low /, = 0.05, ^Icdm = 0.3 , Q\ = 0.65, = 
0, Ho = 65 km/sec/Mpc, Tqmb = 2.726K , Yhc = 0.24, standard recombination). 
The synchrotron spectrum has been calculated assuming for Ei and Bi the scaling 
law (50) with Co = 2.6 (most pessimistic case), and a — 1.5; 1.8; 2.1 and P — 2.9 
respectively. At this frequency (43 GHz) the CMB power becomes comparable 
to the synchrotron power only at very small angular scales (/ > 500). 

Similar calculations at other frequencies confirm that only above ~ 50 GHz 
and at small angular scales (large values of /) the CMBP power spectrum over- 
comes the synchrotron spectrum. Below ~ 50 GHz the CMBP power is always 
below the power level of synchrotron therefore CMBP detection is impossible 
even at small angular scale. 

We can however overcome this limit and plan CMBP observations at lower 
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frequencies if we use our estimator. To see how it improves the CMBP detectabil- 

ity we calculate the angular power spectrum of Df , the estimator we expect when 
the diffuse radiation is completely dominated by synchrotron galactic emission. 
Combining eqs.(39), (44) and (50) we can write: 

This quantity can be directly compared with the power spectrum of CMBP 
because as we saw (see eqs.(31), (37), (40) and (43), Df the power spectrum of 
the estimator evaluated when the sky is dominated by CMBP coincides with the 
CMBP power spectrum. 

Figure 6 shows the power spectra at 37 GHz of \)Df the estimator for a 
galactic synchrotron dominated sky (solid line), calculated using eq.(50), ii)-Df 
the estimator for a CMBP dominated sky which coincides with the power spec- 
trum of CMBP (dotted line) calculated as in figure 5 with CMBFAST using the 
same standard cosmological conditions, iii)the power spectrum of the polarized 
component of the galactic synchrotron radiation (dashed line), calculated using 
eq.(50). As expected at 37 GHz the CMBP power level is well below the level 
of synchrotron power, therefore direct observations of CMBP are impossible (the 
maximum value of the CMBP power is about 2.5 times smaller than the syn- 
chrotron power at the same I). However above I ~ 250 the power associated to 
the synchrotron estimator is definitely below the power associated to the CMBP 
estimator with a maximum ratio CMBP/D ~ 7 at / = 1000. This confirms that 
the use of D allows to recognize the CMBP also at frequencies well below 50 
GHz. 



7 Simulation 

To further test the capabihty of our estimator we studied the separation of CMB 
and galactic synchrotron during observations using, instead of the expected values 

of Di^ as we did in figure 6, measured values simulated by random series of 
numbers, with gaussian distribution ~ A^(0,1), zero mean and unity variance 
(see eqs (54) - (55)). To simulate measurements we generate 2/ + 1 random 
numbers to represent af^ and 21 random numbers to represent af^. We then 
use these values to work out Di using eq.(36). 

Because of their finite angular resolution observations include an average of 
the signal on regions whose angular extension cover a finite multipole interval 
A/. For instance Boomerang data come from regions whose angular extension is 
equivalent to AZ ~ 100 [3], [4], [5]. To get more realistic data we averaged Di on 
an interval /SI — n: 

2 k=+n 

A = rl E A-n/2+.l (52) 
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Figure 6: Power spectrum at 37 GHz of the expected value of the estimator for a 
galactic synchrotron dominated sky (solid line). The dotted line represents both the 
estimator and the CMBP power spectrum for a CMBP dominated sky. The dashed line 
is the power spectrum of the polarized component of the galactic synchrotron radiation 
(dashed line) (see text for details of model and calculations) 
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1 10 100 1000 10000 

Figure 7: Multipole power spectrum of simulated measurements of the estimator D for 
a synchrotron dominated sky (see text) with infinite angular resolution (no smoothing 
on Z, {M = 1)) 

Figure 7, figure 8 and figure 9 give plots of \Di\ (the sign of Di is arbitrary) 
versus /, for A/ = 1, A/ = 10 and A/ = 100 respectively: the very large fluctua- 
tions of the estimator one observes when M — 1, above I ~ 400 are drastically 
reduced as soon as AZ increases. The smoothing effect of the average can be ap- 
preciated comparing figure 7, figure 8 and figure 9 where Di is plotted for n = 1, 
n = 10 and n = 100 respectively: 

Figure 10, figure 11 and figure 12 are similar to figure 10 and plot at 37 GHz, 
instead of the expectation value of Df, simulated measurements of it, for Al — 1 
(figure 10), Al — 10 (figure 11) and Al — 100 (figure 12), respectively. Once 
again the CMBP power spectrum has been calculated with CMBFAST assuming 
the same cosmological conditions we assumed above and the synchrotron power 
spectrum comesfrom eq.(50) with Co = 2.6, a = 1.5,(3 = 2.9 (most pessimistic 
condition) . 

Figure 13 is the same of figure 12 at a much lower frequency, 17 GHz. Here on 
the vertical axis we used logarithmic instead of linear scale, which allows a better 
appreciation of the differences among the three curves. The estimator power 
spectrum now almost touches the two highest peaks of the CMBP spectrum. 
Therefore, 17 GHZ is probably the lowest frequency at which one can hope to 
detect CMBP using Di, in the most favorable conditions. In the most pessimistic 
case (Co = 2.6 and a — 1.5) the corresponding frequency is 25 GHz. 
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Figure 8: Same as figure 7 witli finite angular angular (smoothing on Al = 10) 
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Figure 9: Same as figure 7 with finite angular resolution (smoothing on Al = 100) 



21 



-1 1 



2.0x10 



1.5x10 ' ' - 



o 



CNj 1.0x10" 

+ 



5.0x10 




-12 



1500 



Figure 10: Similar to figure 6. Here we plot simulated measurements with infinite 
angular smoothing resolution (no smoothing on I) instead of the expectation value of 
the estimator at 37 GHz for a synchrotron dominated sky. The dotted line give both 
the estimator and the CMBP power spectrum for a CMBP dominated sky. The dashed 
line is the power spectrum of the galacic synchrotron when its expected contribution 
is maximum (eq. (50) with /? = 2.9, a = 1.5 and Co = 2.6) (see text for details) 
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Figure 11: Same as figure 10 with finite angular resolution (smoothing on AZ = 10) 
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Figure 12: Same as figure 10 witti finite angular resolution (smoothing on Al = 100) 
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Figure 13: Same as figure 12 at 17 GHz using logarithmic scale on the vertical axis 
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8 Conclusions 



Observations of the CMB polarization are hampered by the presence of the galac- 
tic polarized foreground. Only above ~ 50GHz the cosmic signal is definitely 
above the galactic synchrotron and direct observations of CMBP are, in principle, 
possible. Between ~ 30 GHz and ~ 50 GHz the level of the polarized component 
of the synchrotron foreground is at least comparable to the CMBP level, but the 
observational situation is still insufficient to evaluate precisely its contribution to 
maps of the diffuse polarized radiation measured by a telescope. Below ~ 30 GHz 
the galactic signal is definitely dominant. So far a common approach for studies 
of CMBP was to make accurate maps of the diffuse polarized radiation from a 
given region of sky at many frequencies and disentangle the various contributions 
modelling their frequency dependence and spatial distribution. 

We have shown that an alternative way is to take advantage of the differ- 
ent statistical properties of the spatial distribution of the main components of 
the polarized diffuse radiation: CMBP and synchrotron (plus dust) galactic fore- 
ground. By measuring the E and B modes of the polarized radiation we can 
build an estimator which improves the background/foreground ratio by a factor 
sufficient to allow firm recognition and extraction of the CMBP contribution from 
single frequency maps at least down to 25 GHz (17 GHz in the most favorable 
conditions) at angular scales < 0.7° {I > 250). 
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Appendix: stochastic properties of harmonics am- 
plitudes 

Here and overall in paper we suppose that af^ are complex random variables 
which satisfy the probability distribution law: 

(63) 

with variance {\af^\'^) = Ef and {\afj\^) = Bf. 



24 



They have all the propertiers of gaussian variables (below we omitt indexes 
E and B in first and second equations): 

oo 

J p{aim)(Paim = 1 

— oo 

(54) 

oo 

/ aimP{aim)(Paim = 

— oo 
— oo 

(55) 

— oo 

Setting 

oo 

(F) = j Fp{aim)d^aim (56) 

— oo 

with current index E or B, it immediately follows: 

— oo 

(57) 

— oo 
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